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Abstract 



Protein mechanosensitive channels (MS) are activated by tension 
transmitted through the lipid bilayer. We propose a theory of lat- 
eral stress relaxation in a bilayer lipid membrane exposed to ex- 
ternal pressure pulse in the patch-clamp experimental setting. It 
is shown that transfer of lipid molecules into a strained region is 
thermodynamically advantageous due to local decrease of the stress. 
Considered stress relaxation mechanism may explain recent experi- 
mental observations (Davidson and Martinac 2003) of adaptation of 
MscL, bacterial mechanosensitive channel of large conductance, to 
sustained membrane stretch. Lateral stress relaxation in the mono- 
layer, which controls the gating of MscL, triggers thermally activated 
transition of the open channels back to the closed state ("adapta- 
tion"). We evaluate the contribution of the hydrophobic mismatch 
between MS channel and lipid bilayer to the energy barrier separat- 
ing open and closed states. Then, using the MscL thermodynamic 
model (Sukharev et. al, 1999), we estimate characteristic adaptation 
times at room temperature to be of the order of seconds, well in the 
range of the experimental data (Davidson and Martinac 2003). Es- 
timated propagation time of the initial channel-opening stress over 
the whole membrane is 4-5 orders of magnitude shorter. 



I. INTRODUCTION 

The functioning of the protein channels in the cell membranes regulates flows 
of ions in- and out of the cell, thus influencing signal transmission in the neural 
networks (Doyle et al., 1998; Ming Zhou et al. 2001). Mechanical stresses in 
the cell membranes control gating of the protein mechanosensitive (MS) channels 
(Sukharev, 2001), which hence play a role of mechanoreceptors in the cellular or- 
ganisms. Therefore, a theory of the stress propagation and relaxation in the lipid 
bilayer membranes is of substantial interest for the fundamental and practical 
purposes (Yeung and Evans, 1995; Cantor, 1999). One of the biological objects 
convenient for experimental study is the large conductance MS channel (MscL) 
in the inner membrane of the Escherichia coli (E-coli) bacteria. It is possible 
to measure a single channel conductance, an flnd its dependence on the internal 
lateral tension in the membrane. The tension, which causes opening of the chan- 
nel, could be evaluated using video microscopy measurements of the membrane 
curvature formed under the applied external pressure (Sukharev et al., 1999). 

The time-dependence of the membrane's conductance observed in (Hase et al., 
1995) and (Davidson and Martinac, 2003) reveals gradual collapse of the ionic 
currents through the MscL's (called adaptation) within seconds after their open- 
ing. According to existing hypothesis (Sachs and Morris, 1998), this may happen 
due to mutual slide of the lipid monolayers, constituting a bilayer membrane of 
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the &coli. This shde would then cause a relaxation of the channel-opening lateral 
stress inside the membrane. Thus, the slide of the monolayers is induced by the 
same pressure pulse, which leads initially to the opening of the MscL's. 

In this paper we propose a theory of such a stress-relaxation process and 
demonstrate that the slide of the monolayers after application of a constant ex- 
ternal pressure gradient across the membrane is indeed thermodynamically ad- 
vantageous. In Section II we outline the main physical ideas and mechanisms 
considered in this paper. In particular, we apply the theory of the interlayer slide 
dynamics (Yeung and Evans, 1995) in order to evaluate characteristic time of the 
stress-relaxation process. We also show that characteristic time of the MscL adap- 
tation following the stress-relaxation in the membrane is actually much longer. 
In general, this is caused by the necessity of a thermal activation involved in the 
channel adaptation process. Using the MscL thermodynamic model (Sukharcv 
et. al, 1999), we estimate characteristic adaptation times at room temperature 
to be of the order of a few seconds, well in the range of the recently reported 
experimental data (Davidson and Martinac 2003). 

In Section III we introduce a model free energy functional of a bilayer mem- 
brane. We add to the free energy an adhesion energy term, which reflects im- 
portant feature of the patch-clamp experimental setting regarding adhesion of 
the hpid monolayer to the glass wall of the pipette. Then, the self-consistent 
calculation procedure is outlined of the membrane free energy in the semi- and 
complete equilibrium states under a constant pressure gradient in a patch-clamp 
setting. In Section IV we consider results of our numerical calculation of the 
free energy difference between the two successive semi- and complete equilibrium 
conformations of a membrane, i.e. : (1) after application of a constant external 
pressure gradient, but before the interlayer slide; and (2) after the slide of one 
monolayer with respect to the other, leading to the complete equilibrium confor- 
mation. The influence of the adhesion of the lipids to the glass wall of the pipette 
is also studied numerically. Results reasonably compare with the experimental 
data. 

In Section V we review our main estimates of the time scales characterizing the 
membrane and the MscL dynamics. Also, the effect of the energy barrier sepa- 
rating open and closed states of the MscL on its adaptation to sustained stress 
is discussed. The difference between MscL adaptation in the thick and thin lipid 
bilayers is considered in relation with the experimental data for the PC-20 and 
PC-18 liposome patches described in the companion paper (Davidson and Mar- 
tinac 2003). 

Possible experimental verifications and future improvements of the theory are 
also discussed. We focus on the physical ideas in the main text and place more 
detailed mathematical derivations in the appendices. 
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II. BILAYER LIPID MEMBRANE WITH MSCL: CHARACTERISTIC 

TIME SCALES 



A. The shortest time 



The first, semi-equilibrium conformation of a bilayer in the patch-clamp 
setting is achieved shortly after application of a pressure difference pulse P, 
within the time of propagation of the mechanical stress along a membrane: 
At ~ Lo^^ ~ 10~^-=-10~^sec. Here characteristic bending frequency ajj, is evaluated 
below in Eqs. (65) and (64) using the theory of bending waves in a thin plate. In 
the semi-equilibrium state under consideration the number of the phospholipid 
molecules in the curved part of the bilayer membrane has increased relative to 
the initial fiat conformation mainly due to a tearing off the membrane from the 
walls of the pipette, see Fig. 1. In order to allow for this, a finite adhesion energy, 
Ead, is introduced below. The tensions Ti and T2 in the two monolayers of the 
membrane approximately coincide, since the radius R of the membrane is of the 
order of l-=-5/im, while its thickness, 2-Iq, is about (Hasc ct al., 1995) 3.5 ^4nm. 
Hence, on such a short time scale At the membrane behaves effectively like a unit 
structure, as if the monolayers are strongly coupled together and can not shde 
with respect to one another. 



Described conformation is a semi-equilibrium one because it possesses a ten- 
sion gradient along the surface of the lower monolayer 1. The gradient V^Ti 
arises as long as the tension in the curved part of the membrane (in the pipette) 
is finite: Ti 7^ 0, while it is zero in the rest of the patch, see Fig. 1. Hence, layer 
1 will slide against layer 2 (the latter is stuck to the walls). A sliding velocity Vg 
is determined by the drag coefficient (Yeung and Evans, 1995) b : 



where Ka and a are area compressibility modulus for monolayer and dilation 
field respectively (a = a/ao — 1 where a is the area per lipid molecule). Then, a 
conservation equation for the dilation field takes the form of a diffusion equation: 



where D plays the role of a diffusivity constant of the dilation field. Then, the 
lateral stress relaxation time, Tg, which is necessary for the dilation field to level 
off diffusively between the outer patch and the center of the curved membrane 
along the distance ~ is evaluated as: 



B. The lateral stress relaxation time 




at b 



(2) 




i?2 



2R%a 1.2 • lO-^cm^ • Wdyn ■ sec/cm^ ■ 0.05 



= O.Olsec, 



K 




(3) 
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where wc evaluate Ka ^ Ti/a. The characteristic area dilation per lipid, which 
does not destroy the membrane, is evaluated as ci < 0.05. A typical value of h for 
a flat membrane is (Yeung and Evans, 1995) ~ IQ'^dyn- sec/cw?. Measured value 
of external pressure (Sukharev et al., 1999) is about 50mmHg, which corresponds 
to the membrane's lateral tension T — 11. 8dyn/ cm at the curvature radius r ~ 
R ~ 3.5/xm. Hence, the instantaneous tension in the monolayer 1 is: Ti = T/2 ~ 
6dyn/cm, thus leading to above estimate of the lateral stress relaxation time r^. 

The time Tg may become even greater if interdigitation between the monolay- 
ers increases in the curved-stretched conformation of the membrane. The effective 
drag coefficient b may also be enhanced due to an extra impedance caused by a 
finite density of the protein channels piercing the membrane. 

The estimate in Eq. 3 is still considerably shorter than the experimental 
adaptation times ~ Isec (Hase et al., 1995), (Davidson and Martinac 2003). 
In order to understand the reason for such a difference one has to allow for 
additional factor influencing the MscL adaptation time. After the relaxation 
of the lateral tension Ti in the lower monolayer the MscL has to overcome a 
(free)energy barrier to make transition from the open to the closed state. This 
transition, as is estimated below, requires time, which, indeed, proves to be much 
longer than r^. 



C. The MscL adaptation time 

The structure of MscL is assymmetric and the cytoplasmic half of the channel 
poses a barrier to ion permeation (Sukharev et al., 2001). This half of the channel 
is situated in the lower monolayer of the membrane in the considered experimental 
settings (Sukharev et al., 1999; Hase et al., 1995). Hence, activation/closing of 
the MscL is regulated by the lateral tension Ti in this lower monolayer 1. Let us 
use a simple elastic model of MscL channel (Sukharev, Sigurdson et al., 1999) to 
study the kinetics of the channel closing induced by interlayer slide. A detailed 
derivation is presented in the Appendix B (see also Fig. 4). Here we merely 
describe the main consequences. The value of the rate constant for the open to 
closed transition, kc, is the key factor: 

kc^k ex.p{-Eact,c/kBTo} , (4) 

where ks is Boltzmann'constant, Tq is temperature, k is "attempt rate". The 
activation barrier for closing of the channel starting from the open state equals 
(see Appendix B, Eq. (B8)): 

Ea^tAT) -Y^Ao- A,r + T{Ao - A) + ^ (5) 

where T is lateral tension stretching the channel; Aq > A^ designate the areas 
of the pore in the open and top-of the -barrier conformations respectively; Bq 
is elastic constant of the channel in the open state. The parabola Eq. (5) has 
minimum at T < 0. Hence, the activation barrier for the channel closing remains 
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finite (does not reach zero) within the interval of non-negative tensions (T > 0). 
Simultaneously, we find that relaxation of the tension Ti in the lower monolayer 
favours closing of the MscL with the rate kc defined in Eq. (4), as Eact,c in 
(5) decreases with the decrease of T. The pre-exponential factor k in Eq. (5) 
may be estimated as ~ lOOsec"^ (Sukharev and Markin, 2001). Then, regarding 
adaptation time as the equilibration time of the initially open channel Teqv = 
k~^ ~ l.Ssec (Davidson and Martinac, 2003)), triggered by the relaxation of 
lateral stress in the lower monolayer, we estimate the order of magnitude of the 
activation energy Ea,ct,c inverting Eq. (4): 

Eact,c = keTo In k/K = keTo In {kreg^} ^ keTo ln{150} ^ BkeTo . (6) 

We show in Section V that this value of Eact,c has the relevant scale for the known 
set of the thermodynamic parameters characterizing the MscL channels. 



D. The MscL activation time 

The necessity to overcome a free energy barrier may also increase MscL transi- 
tion time from the closed to the open state (MscL activation time) . Remarkably, 
unlike in the case of adaptation, the hight of the barrier for MscL activation 
may vanish at finite tension T* > 0. Thus, the lateral tension T enhances MscL 
opening rate, and the latter process may happen without the thermal activation. 
Using the same elastic model as in Eqs. (4)-(5) above, we find expression for the 
rate constant of the closed to open state transition, ko : 

ko^k exp {-Eact,o/kBTo} , (7) 

where k is "attempt rate" in the closed state. The activation barrier for opening 
of the channel starting from the closed state equals (sec Appendix B, Eq. (B9)): 

Ea^tAT) = f (A - AcY - T{A -^c) + ^, (8) 

where Ab designates the area of the pore in the closed MscL conformation; Be is 
elastic constant of the channel in the closed state. The maximal value Eact,o{T) 
possesses at T = 0: 

E,,,AT = 0) = ^{A-Acr , (9) 

see Fig. 4. The minimal value Eact,o{'^) = is reached at T — T*, where : 

T* = Bc{Ab-Ac); Eact,o{T > T*) = . (10) 

Hence, it follows from Eqs. (7) and (9), (10) that activation time ~ k~^ 
changes from its maximal value Ta{T = 0) at zero tension down to its shortest 
value at T > T*: 
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k-'e^p{Bo{Ab-Acf 
k-\T> T* 



/2kBTo} , T = 



Indeed, according to (Sukharev, Sigurdson et al., 1999), an increase of T from 
11.9mN/m to lA.lAmN/m leads to a decrease of the activation time Tq ~ by 
about 40 times, so that it reaches the absolute value of 6.8 • 10~^sec. 



Wc have also calculated the free energy of the whole patch: bent part and 
unstrained outer part, as a function of the external pressure difference P in both 
the semi- and complete equilibrium states. Indeed, the free energy in the relaxed, 
Ti = state, proves to be lower than in the semi-equilibrium state achieved before 
the slide. Fig. 2 a. A comparison of the calculated numbers Ni and N2 of the 
lipids, in the curved lower and upper monolayers correspondingly, in the semi- 
and complete equilibrium states follows from Fig. 2b. It shows that the number 
ANi of the lipids "sucked in" from the patch to relax Ti is of a macroscopic 
magnitude, e.g. AA^i ~ 4:%Ni. Based on these result, we estimate also a diffusion 
time, Tfi, which would be necessary for the lipid molecules to move from the patch 
to the monolayer 1 chaotically, i.e. without the slide motion of the monolayer as 
a whole. Experimentally determined diffusion coefficient of the individual lipid 
molecules in bio-membranes at room temperature is (Sonnleitner et al., 1999): 
D ~ 1 lOiim^ / sec. Hence, we evaluate a characteristic diffusion time as: 



which is by orders of magnitude longer than the sliding time r^, Eq. 3. Hence, 
diffusion of individual lipid molecules proves to be less effective for the relaxation 
of the stress in the concave monolayer than the interlayer slide. 



Despite the seeming clarity of the mechanisms, which determine the different 
time scales described above, the whole picture of MscL response to the sustained 
lateral tension in the membrane can be more involved: hydrophobic mismatch 
affects MscL gating (Perozo et al., 2002). While MscL exhibited adaptation in 
bilayer formed with diecosenoyl phosphatidy choline, PC-20 molecules, no adap- 
tation to sustained stress was observed in bilayer made of dioleoyl phosphatidy- 
choline, PC- 18 (Davidson and Martinac 2003). Below we present some estimates 
demonstrating that the change of MscL hydrophobic mismatch between its closed 
and open conformations can be responsible for a fairly different gating behavior 
depending on whether the lipid bilayer is "thick" or "thin". This may explain 
experimentally observed drastic difference betwee MscL gating behaviors in PC- 
20 and PC- 18 bilayers. We use here again the elastic model of MscL channel 
(Sukharev and Mrkin, 2001) in combination with the lipid-protein hydrophobic 



E. The lipid diffusion time 



Td - {ANi/Ni)R^/D - 0.1 Isec, 



(11) 



F. Adaptation in PC-20 versus no adaptation in PC-18 
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mismatch model (Fournier, 1999) shortly introduced in the Appendix. 

The length of hydrophobic region of (the external wall) MscL channel is differ- 
ent for its closed and open conformations (Sukharev et al., 2001). The mismatch 
energies in the open and closed states of the MscL channel, F£p and, F£p, con- 
tribute to the total free energy of these conformations, Fq and Fc- The free 
energy difference between open and closed states is then given by: 

Fo-Fc^Eo + AFlp , (12) 

here AFlp = F£p — F^^p is the part of the energy difference arising from the 
difference in hydrophobic mismatch energies, while Eq is determined by other 
factors (Sukharev, Sigurdson et al., 1999), e.g. the type of lipids in the bilayer 
(hydrocarbon chains lengh; curvature stress profile), thermodynamic conditions 
of the surrounding solution (pH, temperature). In the elastic model of MscL at 
T = 0: Eq + AFlp = Eq. The mismatch energy depends on the type of lipid; if 
the length of the lipid tail provides better matching of lipid-protein hydrophobic 
regions, this results in lower mismatch energy and lower total energy of the given 
channel conformation. In other words, lipids surrounding the channel favour 
certain (closed or open) conformation of the MscL (Hamill and Martinac, 2001). 

For example, if in the (initial) closed state of MscL the hydrophobic thickness 
of the bilayer is equal or greater then that of the channel, the channel opening in- 
creases the mismatch and, thus, the mismatch energy adds to the value Eq (and to 
the energy separation Fq — Fc). On the other hand, if the hydrophobic thickness 
of the bilayer is less or matches that of MscL in the open conformation, the chan- 
nel opening decreases the mismatch and, hence, the mismatch energy subtracts 
from Eq. In the thermodynamic model of MscL it corresponds, see Fig. 4, to the 
decreased vertical distance (i.e. £'o at T = 0) between the parabolas describing 
the free energies of the open and closed conformations of the channel. When 
MscL opens, its hydrophobic thickness decreases by 6 -h 8A (S. Sukharev, Mary- 
land University, personal communication, 2003). Our estimate (see Appendix B) 
shows, that the difference of mismatch energies between open and closed states 
of MscL, AFlp, may range up to 20kBT. 

We would expect that shorter lipids (e.g. PC-18 compared to PC-20) favour 
open state of MscL: 

(AF^p)^^i« < (AF^p)^^2° . (13) 
This implies that the energy difference, Fq — Fc, becomes smaller: 

{Fo - Fcf'^'' < {Fo - Fcf'^''^ . (14) 

Thus at any tension the ratio of the equilibrium fractions of the open and closed 
channels, no/nc ■ 

no/nc ~ exp {-(Fq - Fc)/kB%} (15) 
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is shifted towards increased fraction of open channels no for PC-18 as compared 
to PC-20 : 



(no/nc)^^^^ > (no/nc)^^^° . (16) 

This is why in PC-18 hposome an interlayer shde (and resulting relaxation of the 
lateral tension in the lower monolayer) may not lead to a complete closing back 
of the MscL's in the final thermodynamic equilibrium state, since at Fq ~ Fc 
we have no/nc ~ 1- We also expect the closing of (a fraction of ) MscL's to 
take longer time in PC-18 bilayer in comparison with PC-20 bilayer as the energy 
barrier between closed and open conformations, Eact,c, increases with a decrease 
of Eq, Fig. 5 (see explanation in Appendix B after Eq. (B9)). The higher energy 
barrier leads to the lower rate of closed-to-open transition, kc in Eq.(4). 



G. Interlayer slide and hydrophobic mismatch 

Interlayer slide makes no significant contribution to the energy of hydrophobic 
mismatch. Due to the incompressibility condition,!^ — a ■ I — const {V - is 
volume of lipid molecule, a - is area per lipid molecule, l-is lipid length in the 
strained membrane), tension in each monolayer, Tj, produces fractional change 
in its thickness: 

T,^-^.Ka. (17) 

Here i = 1, 2 is monolayer label; A/ = / — /q, and Iq is the length of a lipid in the 
unstressed fiat membrane; Ka is the area compressibility (expansion) modulus 
of a monolayer. After the application of external pressure, P, to the membrane 
(but before interlayer slide) both monolayers are strained and tensions in them 
approximately coincide: 

Ti ^2 ~ T/2 , (18) 

here T is tension in the membrane defined according to Laplace law: T = Pr/2, 
r is curvature radius. The total relative change in membrane thickness, is given 
by: 

Ml + M2 



2ln 



2Ka 



+ 



T2 
2Ka 



T 

(19) 



2Ka 



After the interlayer slide the flow of lipids from the unstrained outer patch to 
the lower monolayer leads to the relaxation of lateral stress in it; tensions are 
redistributed: 

ri = o,r2 = r. (20) 

But the total change in the bilayer thickness, evidently, remains nearly the same: 

A/2 T 

A 1 = 0, — -p^ . 21 

' 2 2Ka ^ ' 
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III. THE MODEL FREE ENERGY OF MEMBRANE 



We start from the microscopic model of a bilayer membrane as described e.g. 
in (Safran, 1994; Ben-Shaul, 1995; Xiang and Anderson, 1994). First, write a 
single i-th monolayer (i = 1,2) free energy per lipid molecule ff. 

fi = fsi + //^^ + fti = lai + C/ai + {ks/2){k - Isf (22) 

Here fsi is the external surface energy of the monolayer, which increases with the 
area per lipid molecule due to the hydrocarbon tails interaction energy with 
the solvent. The short-range repulsion between the polar molecular heads at the 
surface of the membrane is represented by the second term in Eq. 22. The third 
term, though looks like a spring elastic energy (with not stretched spring length 
Is), codes for the entropic repulsion between the hydrocarbon tails in the depth of 
each monolayer (Safran, 1994). In what follows we omit the second term in Eq. 
22 on the empirical grounds (Ben-Shaul, 1995) relevant for the chain molecules 
in phospholipid bilayers. There is no interface energy term included in Eq. 22 for 
the internal surfaces of the monolayers forming the bilayer, as the latter are not 
reachable for the solvent molecules. We also assume a vanishing interdigitation 
of the tails between the adjoint monolayers. The length of a lipid molecule /j, 
i.e. the thickness of the i-th monolayer, is not an independent variable from the 
aria per molecule a^. The volume v per lipid molecule is conserved, i.e. the latter 
could be considered as incompressible (Safran, 1994; Ben-Shaul, 1995). Due to 
this conservation condition the length is related to the mean and gaussian 
curvatures H and K at the interface between the monolayers by a well known 
differential geometry formulas (Safran, 1994) : 

/i = ^oi-ioi^+(2/3)4X (23) 
l2^k2 + ll2H+{'2/?,)ll^K (24) 

Here the change of sign in front of iJ-term in the equation Eq. 24 relative to 
that in equation Eq. 23 is due to a simple geometrical fact that the external 
normal vector to e.g. the layer 1 is simultaneously an internal normal vector 
to the layer 2 at their mutual interface. Different Zqi parameters just reflect the 
incompressibility of the lipid molecules mentioned above: 

loi = v/ai (25) 

Substituting Eqs. 23 and 24 into the third term of Eq. 22 one finds contributions 
to the free energy from the tails in the form: 

fti = {hltJ2) [{H - Coi)^ + (4/3)CoiZoii^] (26) 

ft2 = {kslo2/'2) [{H + C02)' + (4/3)Co2/o2i^] (27) 

where parameters Coi have the meaning of the local spontaneous curvatures of 
the i-th layer: 
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C'oi — {loi — h)/^Oi 



(28) 



We readily recognize the Helfrisch's formula (Zhong-Can and Helfrich, 1987; 
Zhong-Can and Helfrich, 1989) for the free energy of a membrane in Eqs. 26 and 
27. 

In a symmetric bilayer case the numbers of lipids in the 1-st and 2-nd mono- 
layer are equal: A^i = A^2- Under the equivalent conditions on the opposite 
surfaces of the membrane: Oi = a2, the linear in H terms in Eqs. 26 and 27 
would cancel in the total free energy F: 

^ = E = E Nifi (29) 

1=1,2 i=l,2 

Nevertheless, as shown below, a linear in H term may arise when the up vs 
down symmetry of a membrane is broken, e.g. by applied pressure gradient in 
combination with the different boundary conditions at the monolayers peripheries, 
see Fig. 1. 

To make the whole idea transparent we restrict our present derivation to the 
case of a "spherically homogeneous distributions" of molecules, i.e. considering 
a^'s as being different for the different indices i's, but position-independent within 
a curved part of the i-th monolayer. Below, we neglect the inhomogeneous dis- 
tribution of the strain across the thickness of each monolayer. Also, we consider 
only spherical shapes of the curved part of the bilayer membrane, while intro- 
ducing position independent (over the membrane's surface) mean and gaussian 
curvatures: 

H=l/r; K = H^ = l/r^ (30) 

where r is the radius of curvature. The radius of the base of the curved membrane 
is fixed at i? < r, in accord with the fixed radius of the pipette, which sucks in 
the membrane, and creates a pressure difference P between inside and outside 
surfaces of the membrane in the experimental setup (Sukharev et al., 1999; Hase 
et al., 1995). Then, according to the Laplace's law we have: 

Ti + T2 = Pr/2 , (31) 

where T; is the lateral tension in the i-th monolayer, and we neglected a small 
difference between the curvatures of the monolayers. 

First, we consider the free energy of a flat, undeformed bilayer membrane, in 
order to use it as a reference value of the free energy per lipid molecule, /q. The 
flat membrane has a monolayer thickness Zq ~ v/ao, where the area per molecule 
flg is determined from the condition of the minimum of the free energy F deflned 
in Eq. 29: 

dF/dao = , at: (32) 
01 = 02 = 00, H^K = 0, Ni = N2^No. (33) 

Using conditions from Eq. 33 in the Eqs. 22 - 29 we derive the following expression 
for the free energy per lipid molecule, /o, in the undeformed state: 
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/o = fi{zo, /i = 0) = fiizo, /i = 0) = esZo + et{l + Zq'^ - 2zq ^} , 



(34) 



where the following dimensionless parameters have been introduced: 

Zi^Qils/v; h^ls/r = Hls] ho^ljR; es^^v/lg] et^ksll/2. (35) 

It is useful for the comprehension of the rest of the paper to mention here orders of 
magnitude of the main parameters. Using the values of R and Ig ~ Iq mentioned 
in the Introduction, we find that the dimensionless curvature is of the order: 
h ~ 10~^ <^ 1. Therefore, it could be inferred from Eqs. 23 and 24, and 
from the definition Eq. 25, that v ^ Istti, which leads to an estimate: Zi ~ 1. 
The entropic nature of the fa term in Eq. 22 is reflected in the temperature 
dependent coefficient kg, which for the known phospholipid bilayers is of the 
order (Ben-Shaul, 1995; Xiang and Anderson, 1994): kj^ ~ 25kT, where k is 
the Boltzman's constant and T is the absolute temperature. The scale of the 
surface term fsi in Eq. 22 is characterized by the coefficient 7, which at room 
temperature is of the order (Ben-Shaul, 1995): 7 ~ O.lkT/A^. Gathering all 
the estimates, and also taking into account that the typical value of an area per 
molecule is (Ben-Shaul, 1995): Oj ~ 60A^, we obtain the following list of the 
estimates (taken for the room temperature T ~ 300K) : 

Zi^l; 10"^; ~ 3 • lO^^Vg; ~ 1.3e,. (36) 

An additional important dimensionless parameter relates characteristic experi- 
mental value of the pressure difference P ~ 50mm Hg (Sukharev et al., 1999; 
Hase et al., 1995) with the "microscopical bending energy" et- 

Pv/{2et) ~ 10"^ at: ~ 5 • lO'^^erg; -y ~ 1.2 • IQ^A^, (37) 

where we use an estimate Ig ~ 20A for the typical length of a free phosphohpid 
molecule (Ben-Shaul, 1995) and = O.Gej. 

In the dimensionless parameters and variables equation Eq. 32 reads: 

dh/dzo = £3 + 2et{z^^ - z^^] = , (38) 

or in the Cardano's form: 

zl + PqZq + g = 0; where: pQ = -q = 2et/es (39) 
and: {q/2f + (po/3)' = (e^/e,)' + (2ei/3e,)' > 0. (40) 

The inequality in Eq. 40 proves that there is a unique real root of the cubic 
equation Eq. 39. Solution zq depends only on the dimesionless parameter po = 
2tt/ts. A numerical solution of Eq. 39 at 1/po = 0.3 gives: 

zo = 0.829 . (41) 

Substitution of the solution zq into Eq. 34 gives desired value of /q. 
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A. Free energy including patch and adhesion to the walls 



The free energy expressions Eq. 22 and Eq. 29 do not account neither for 
the adhesion of some part of the upper monolayer to the wall of the pipette, nor 
for the existence of the patch, which serves as a reservoir of lipids for the lower 
monolayer, see Fig. 1. In order to include these features of the experimental 
setting (Sukharev et al., 1999; Hase et al., 1995) we write the total Helmholtz 
free energy of the membrane as follows: 



F = iVi(/i - /o) + N2{h -fo + Ead) - iVi(ai - aoi)T, - N2{a2 - 002)^2 - 



Explain the right hand side (rhs) expression in Eq. 42 term by term. The first 
two terms signify the free energies of lipids in the curved parts of the monolayers 
1 and 2 correspondingly. The energy per lipid molecule /o is subtracted in order 
to allow for the free energy difference between the molecules which belong to the 
curved part and to either the patch, or to the part of the monolayer 2 that is stuck 
to the wall of the pipette. In the latter case, we allow also for the adhesion, so that 
the energy per molecule is lowered with respect to the patch by an amount of the 
adhesion energy Ead- The free energies per lipid molecule fi in the i-th monolayer 
could be written using definitions Eqs. 22-29 and dimensionless variables Eq. 35 
as follows: 



h) = e,^i + et {(/iV3)(7^r^ - 4z^') + 2h{z^' - z^') + {z^' - 1)'} ; (43) 



The next two terms on the rhs of Eq. 42 represent a mechanical work done 
by the tensions Tj=i^2 by stretching an area per molecule in the curved part of 
the i-th monolayer from the initial to the final value, and a, respectively. The 
initial areas Ooi, corresponding to the equilibrium state at zero external pressure 
gradient, are calculated below together with the initial numbers A^oi of the lipids 
in the curved parts of the monolayers. The fifth term on the rhs of Eq. 42 is 
introduced to allow for the mechanical work done by the perpendicular to the 
wall projections of the tensions during tearing off the 2nd monolayer from the 
pipette's wall. Angle © is the wetting angle between the membrane and the wall: 



sin e = {1 - cos' e}^/2 = {1 - {R/rfy^^ = {1 - {h/hofY'^ , (45) 



where notation defined in Eq. 35 is used. The factor sin © • (Ti + T2) in the fifth 
term signifies that the force acting on the 2nd monolayer perpendicular to the 
wall consists of the sin © ■T2 projection added to the projection of sin © -Ti, which 
is applied from the side of the 1st monolayer. The latter force equals, with an 
opposite sign, to the force acting perpendicular to the wall on the monolayer 1 in 
accord with the Newton's 3rd law. Next, the terms with the Lagrange multipliers 
\=i,2 are introduced in order to select only the spherically shaped conformations 
of the membrane as its curvature H changes together with the external pressure. 



(A^2 - A^o2)aosin© • (Ti + ^2) - E USi{H) - N,ai). 



(42) 



f2{z2,h) = fi{z2, -h). 



(44) 
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In this case the external surface areas Si and S2 of the monolayers and the 
numbers of the molecules in them are related by the following equations: 

S,{H) = 27r(r - hf (l - {1 - R'/ir - hYy/') = N,ai ; (46) 

S2iH) = 27r(r + (l - {1 - i?V(^ + kfV^^) = N^a^ . (47) 

Here the thicknesses /j=i,2 of the i-th monolayer are defined in Eqs. 23 and 24. 
As long as the ratios are small: 

k/r<k/Rr^lO-\ (48) 

we shall neglect small corrections k to the curvature radius r in the expressions 
in Eqs. 46 and 47, hence defining a unique surface area function S{H): 

S{H) = 2nr^ (l - {1 - R^rY^^) = 2t:H-^ (l - {1 - {HRfY^^) (49) 

Then, using Eq. 49 and relations Eq. 46, Eq. 47, together with definitions Eq. 
35 we express the number of the lipids in the fiat bilayer Nq as follows: 

So = S{H ^ 0) = ttR^ = Noao = NqZqv/1, . (50) 



B. Equilibrium conformation under zero external pressure difference: P = 

In the absence of pressure, the membrane is essentially flat, subjected to the 
"resting" (Sukharev et al., 2001) tension arising from the membrane adhesion to 
the (glassy) surface of the pipette. Thus, consider first a bilayer membrane in the 
equilibrium state in the pipette before an external pressure gradient is applied. 
Then, the tensions induced according to the Laplace's law Eq. 31 by the external 
pressure difference P are zero: 

Ti = = , (51) 

and the free energy equals: 

F = iVoi(/i - /o) + A^02(/2 - /o + Ead) - E HS{H) - NoiUoi). (52) 

1=1,2 

Corresponding equilibrium equations minimizing F in Eq. (52) are derived in 
Appendix. Using also relations Eq. 46 and Eq. 47, we obtain the following 
closed form equation for the unknown value of H expressed in the dimensionless 
variable h defined in Eq. 35: 

E nM/dh = -{l/No){dS{h)/dh) dfi/dzi ; (53) 

i=l,2 i=l,2 

where: m = N,/No = 2{zo/ z{){ho/hf (l - {1 - {h/hofY'^) . (54) 
The functions fi are defined in Eqs. 43 and 44. 
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C. Semi-equilibrium conformation at finite external pressure difference P 

Soon after an application of a finite external pressure difference P perpen- 
dicular to the plane of the membrane, i.e. during the times < i -C r^, a 
semi-equilibrium conformation is achieved. The free energy of the membrane in 
the semi-equilibrium state could be written using Eq. 42 as follows: 

F^F- Xs{Ni -N2- iVoi + N02) . (55) 

Hence: 

F = m{f, - /o) + iV2(/2 - /o + Ead) - (iV2 - iVo2)ao sin Q ■ {T, + T^) - 
J2 [Niiai - aoi)T, + K{Si{H) - TV-a,)] - \-,{N, - N2 - N^i + N02) . (56) 

i=l,2 

Here an extra Lagrange multiplier A3 is introduced to reflect a simple physics. 
Namely, as explained in the Introduction, we consider the semi-equilibrium state 
as existing during such a short time, which is not enough for a transfer of the 
lipid molecules from the patch to the curved part of the monolayer 1. Hence, a 
change of the numbers of lipids Ni{P)—NiQ, in the curved parts of the monolayers 
i — 1,2 relative to the initial, P — equilibrium conformation occurs only due 
to a tearing off the 2nd monolayer from the pipette's wall. In the latter case it 
is reasonable to assume that the numbers of lipids in both monolayers increase 
simultaneously by the same amount (no interlayer slide): 

N, - Noi ^N2- N02 . (57) 

Here A^oi are involved, which had to be found from the Eqs. A1-A4. Condition Eq. 
57 is maintained by the term with A3 in Eq. 56. As explained in the Introduction, 
due to small thickness of the membrane Eq. 48 we assume equal lateral tensions 
in the monolayers, which then could be found from the Laplace relation Eq. 31: 

Ti = T2 = Pr/4 = T/2 , (58) 

Here new parameter T is introduced for convenience. Corresponding equations 
for the minimum of Eq. 56 are presented in the Appendix. 



D. Relaxed conformation at finite external pressure difference P 

Equations for a complete (relaxed) equilibrium conformation of the membrane 
are readily obtained, see Appendix, starting from the following expression for the 
free energy: 

F = A^i(/i - /o) + A^2(/2 - /o + Ead) - {N2 - No2)aoTsmQ - 

N2{a2 - ao2)T - ^ Xi{Si{H) - N^a,) . (59) 

i=l,2 
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Equation 59 is obtained from Eq. 56 by "relaxing" to zero the tension Ti in the 
monolayer 1, and by dismissing a condition of a "dynamic cut off' of the first 
monolayer from the reservoir of lipid molecules, i.e. from the patch. The latter 
condition was represented in Eq. 56 by the term with the Lagrange multiplier 
A3. The lateral tension in the monolayer 2 obeys the following Laplace relation: 

T2 = Pr/2 = T ;Ti = 0. (60) 

Corresponding equations for the minimum of Eq. 59 are presented in the Ap- 
pendix. 

IV. DISCUSSION OF NUMERICAL RESULTS 
A. Initial equilibrium state 

Results obtained by the numerical solution of equations A1-A4 are represented 
in Table 1. Despite the absence of the external pressure gradient: P = 0, the 
curvature 1/r proves to be small but finite: R/r ~ 10~^. This happens due to 
finite adhesion energy Ead of the monolayer 2 to the wall of the pipette. In Table 
1 n2 is the number of molecules in the curved part of the monolayer 2 normalized 
with the number Nq in a free monolayer. Presented data indicate that 712 decreases 
with an increase of the adhesion energy Ead- The relative decrease of n2 ranges 
from 0.8% at Ead/^t = 0.03 to 9% at Ead/^t = 0.3. Simultaneously, the relative 
number of molecules rii (ni = Ni/Nq) in the monolayer 1 remains practically 
constant and equals 1. The reason for the different behaviour of rii and n2 is 
obvious. More molecules from the curved part of the monolayer 2 tend to stick 
to the wall as the adhesion energy Ead increases. On the other hand, though 
the monolayer 1 leans to the wall following the shape of the monolayer 2, it does 
not stick, hence it compensates for the "losses" of the number of molecules in its 
curved part by absorbing extra molecules from the patch. 

Table 1 also contains the values of the areas ai^2 per molecule in the mono- 
layers, expressed in dimensionless units via the variables zi, Z2 defined in Eq. 35. 
The range of Ead/^t is chosen from the consideration of the stability of the mem- 
brane under stretching of its area, which is known to have an upper threshold 
(Sukharev et al., 1999) of about (oj — ao)/ao ~ 4%. Remarkable is that in this 
way our theory predicts reasonably well the value of the adhesion line tension 
of a lipid membrane to the glass wall known from the experiment (Opsahl and 
Webb, 1994). Indeed, theoretically chosen ratio Ead/^t = 0.03 in combination 
with experimentally known absolute values ej ~ 4 • 10~^^ergi and oq ~ 60^4 leads 
to the estimate : 

Ta ~ Ead/aQ = 2.5 dyn/cm , (61) 

which falls right inside the experimentally measured interval Ta — 0.5-^ A dyn/cm 
(Opsahl and Webb, 1994). 
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Finally, the last line of Tabic 1 provides normalized values of the free energy of 
the membrane F, defined in Eq. 52. It is interesting to notice that the values of 
F/ (iVoej ) practically coincide with the corresponding values of the (normalized) 
adhesion energy per molecule E^dj^u which are also indicated in Table 1. An 
examination of Eq. 52 leads to the conclusion that the coincidence is not acci- 
dental. Namely, while the reference of energy for the molecules in the monolayer 

1 is approximately /o per molecule, the one for the monolayer 2 is /o — Ead- At 
zero external pressure P = the free energy per molecule in each monolayer ap- 
proximately equals /o, i.e. the fiat membrane's value. Therefore, contribution of 
the monolayer 1 to F in Eq. 52 is nearly zero, and contribution of the monolayer 

2 is indeed EadNg, as long as A^o2 ~ A^o- 



B. Conformations at P 7^ 

Numerical results were obtained by minimizing the free energy expressions 
of the memebrane in the pipette under a fixed external pressure. Expressions 
Eq. 56 and Eq. 59 were used, describing correspondingly the semi- and complete 
equihbrium states of the membrane. The calculated theoretical dependences are 
presented in Figs. 2, and 3. Based on the experimental data (Sukharev et al., 
1999; Hase et al., 1995; Sukharev et al., 2001) we had chosen the following values 
for the two dimensionless ratios controlling the main dependences: 

e,/et = 0.6; F^^e* = 0.03 . (62) 

The first ratio in Eq. 62 sets the scale of the surface tension, while the second 
ratio fixes the scale of the adhesion energy of the membrane to the glassy wall 
of the pipette with respect to the energy of the entropic repulsion between the 
hydrocarbon tails in the depth of each monolayer. 

As is shown in Fig. 2 a, the free energy of the relaxed conformation is lower 
than that one in the semi-equilibrium state at any pressure difference P measured 
in the dimensionless units introduced in Eq. 37: 

P ^p = Pv/{2et). (63) 

The Fig. 2a justifies our conjecture made in the Introduction that a relaxation to 
zero of the lateral tension Ti in the concave monolayer 1 is concomitant with the 
decrease of the free energy F of the bilayer membrane kept at a fixed external 
pressure difference. A slide of the monolayer 2 is prevented by its finite adhesion 
enrgy to the wall of the pipette, i.e. Ead- 

Next, Fig. 2b provides an important evidence that the decrease of the free 
energy and vanishing of Ti happens by virtue of an increase of the number of 
lipid molecules in the curved part of the monolayer 1 at a constant pressure 
difference P. Indeed, n^""" is manifestly greater than nf'^'^""' at the same pressure. 
Simultaneously, nf'"" is even shghtly smaller than nf-^'^. The increase of n2 under 
fast growth of the external pressure P happens due to a tearing off the membrane 
from the glassy wall of the pipette. In the absence of the transfer of lipids from the 
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patch, an "instant" increase of n^'^ with pressure in the semi-equihbrium state 
foUows that one of nf"'^'^ and has the same origin, i.e. tearing off the membrane 
from the wall. But, after enough time had elapsed since the beginning of the 
pressure pulse, the rii shows an increase from nf^'^ up to at a fixed pressure 
due to a transfer of lipids from the patch. The latter effect happens either via 
interlayer slide or/and via lipid diffusion, as discussed in the Introduction. During 
this relaxation process the area per molecule in the monolayer 1 decreases from 
^disc j^jj^ dimensionless units) to practically that one in the free bilayer: zf""'"' ~ 
Zq, see Fig. 3a. Simultaneously, the monolayer 2 has to bear now the whole lateral 
tension: T2 — T; Ti — (previously shared with the monolayer 1: T2 — Ti — 
T/2). The area per molecule in the monolayer 2 stretches from z***^ to ^2"""- As 
long as an increase of this area relative to a free bilayer value exceeds 4 ^ 5% 
the membrane becomes unstable. Therefore, we deduce from Fig. 3a, that the 
range of the external pressure differences in our model should be bound as follows: 
< P < 1.5. 

Results of our calculations of the membrane's curvature as a function of the 
pressure P are presented in Fig. 2 c. There is apparent effective softening of the 
membrane after a relaxation at constant P, i.e. h,.on,i > hdisc- Also, the curves 
exhibit tendency to a saturation of the curvature at hight enough pressures, in a 
qualitative accord with the experimental data (Sukharev et al., 1999). Neverthe- 
less, according to our results plotted in Fig. 3b, the lateral tension Ti continues to 
grow smoothly as a function of pressure P even in the region of (near) saturation 
of the curvature. 

Finally, we mention here, that separate consideration of the semi- and com- 
plete equilibrium states discussed above, is justified by the fact that the charac- 
teristic stress relaxation time: Tg^d ~ 0.1-=-lsec differs by 4-=-5 orders of magnitude 
from the characteristic time At necessary for a propagation of the bending defor- 
mations over a membrane with the linear dimensions R ~ Ifim: 

At ~ {ph/etf^ ~ 10~^cm^ {ig/cm^ • cm/lQ-^^ergf'^ ~ IQ-^sec . (64) 

Here we used a well known formula for the frequency of the bending waves in a 
thin plate (Landau and Lifshitz, 1980): 

uJb = q^{kB/{2lop)Y'\ (65) 

where /c^ is the bending modulus, g is a wave- vector, p is the mass density of the 
membrane, and 21q is its thickness. Actually, the latter is dynamically increased 
by a drag of the adjacent layers of the liquid in which the membrane is immersed. 

In Eq. 64 we use as an estimate of fcg involved in Eq. 65 allowing for the role 
played by ClS Si bending coefficient. This is justified by a direct comparison of 
the free energy expressions Eqs. 26, 27 and the definition of et given in Eq. 35. 

V. SUMMARY 

Dynamics and characteristic times of bilayer membrane with MscL channels 
in the patch-clamp experimental setting are studied theoretically. The different 
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characteristic times range from 10~^ -7- 10~^ sec to several seconds. The shortest 
time corresponds to stress propagation along the micrometer-sized membrane 
bent by external pressure gradient. The external pressure produces lateral tension 
in the bilayer, which opens mechanosensitive channels (MscL). We show that 
lateral tension in lipid bilayer does not remain constant under constant pressure 
gradient applied to the membrane. Interlayer slide leads to redistribution of 
lateral stress within 0.01^0. Isec after application of the pressure step. Relaxation 
of lateral stress in the lower monolayer triggers closing back of MscL channels 
("adaptation"), observed experimentally, see companion paper. The longest time 
(~ 1 sec) arises from the free energy barrier for the transition from the open to 
the closed MscL conformation. 

We have analysed the influence of bilayer hydrophobic thickness on the hy- 
drophobic mismatch energy and MscL thermodynamics. Thick bilayer hinders 
MscL opening, while thin bilayer favours the open conformation of the channels. 
Our estimates of the mismatch energy give values in the range 2 20A;sTo, where 
ksTo is the thermal energy at room temperature Tq. This estimate is compa- 
rable with the free energies in the MscL elastic model (Sukharev et al., 1999) 
and (Sukharev and Markin, 2001). Depending on the bilayer thickness, the hy- 
drophobic mismatch favours either closed or open MscL conformation, as outlined 
in Section IL This may explain why MscL adaptation phenomenon, being present 
in the PC-20 lipid bilayers, is absent in the PC-18 bilayers (see companion paper). 

We have studied the free energy functional of the membrane modeling the 
patch-clamp experimental setting. Minimizing the free energy functional, we have 
found parameters of the membrane conformation as function of pressure gradient 
for two (semi- and complete) equilibrium states. The semi-equilibrium state is an 
instantaneous mechanical equilibrium state achieved after application of external 
pressure step to the membrane. Following it, a transfer of lipid molecules along 
the lower monolayer (see Fig. 1) takes place. It leads to the leveling of the 
finite lateral tension in the strained part of the lower monolayer with zero tension 
existing in the outer patch. As a result the free energy reaches its true minimum 
and the membrane achieves the complete equilibrium state. The upper monolayer 
is fixed due to lipid adhesion to the pipette glass wall. Our free energy functional 
of the membrane includes this adhesion energy. Theoretically predicted value of 
the adhesion line tension in our model is in a good agreement with experimental 
data (Opsahl and Webb, 1994). 

The transfer of lipids in the lower monolayer may occur either via interlayer 
slide or/and via diffusion of lipids into the stretched region of the lower monolayer 
from the outer patch (see Fig. 1), see also (Baoukina and Mukhin, 2003) and 
(Mukhin and Baoukina, 2002. Lipids diffusion mechanism of stress relaxation 
in a bilayer fluid membrane under pressure, cond-mat/0206099). Characteristic 
time of the individual lipid diffusion processes is by one order of magnitude longer 
than the time of the interlayer slide (collective motion of the lipids constituting 
the lower monolayer). Nevertheless, this times may become comparable at high 
enough density of MscL's, which then would pin the monolayer against slide. 

The above explanation of the slide-triggered closing of a MscL leads to a 
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proposal (Th. Schmidt, Leiden University, personal communication, 2002) of a 
possible experimental check up of the theory. Namely, suppose that the membrane 
with MscL has been turned over, so that the cytoplasmic half of the channel 
belongs now to the monolayer 2, which in turn stucks to the wall of the pipette. 
Then, the opening of the channel would be controlled by the lateral tension T2, 
which does not vanish (it even increases) after the relaxation of the membrane 
under a constant pressure. Hence, one would expect that adaptation of MscL 
after an application of the opening pressure would be less probable in this case. 

Finally, we mention some possible future improvements of the theory. These 
could be made by lifting the simplifying restrictions of the homogeneity of the 
lateral stresses (tensions) across the thickness of each monolayer. Besides, one 
may abandon the approximation of a position-independent area per molecule 
on the membrane's surfaces. As a consequence, in a more elaborate scheme 
of derivations one would not restrict himself to the spherical symmetry of the 
membrane's shape. Nevertheless, we believe that these improvements would not 
deny the main physical concepts considered in this paper. 

The authors are grateful to S.L Sukharev for the introduction in the problem and 
for numerous enlightening discussions during the work. Useful discussions with Jan 
Zaanen and comments by R.F. Bruinsma and Th. Schmidt are highly acknowledged. 
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APPENDIX A: EQUILIBRIUM EQUATIONS 
1. Equilibrium conformation, P = 



Considering Nq^, Ooi and H as the five independent variables, we minimize F 
in Eq. (52) with respect to all of them, and thus obtain the following equilibrium 
equations: 

dF/dNoi = 0: {/i - /o}/aoi = -Ai ; (Al) 
dF/dNo2 = 0: {f2-fo + Ead}/ao2^-X2; (A2) 
dF/daoi = 0: dfi/daoi = -Xi ,i = 1,2; (A3) 
dF/dH = : Yl Noidfi/dH = (Ai + \2)dS{H)/dH . (A4) 

1=1,2 

One can easily see that Eqs. Al and A2 permit exclusion of the Lagrange mul- 
tipliers Aj from the equations A3 and A4. Then, solving Eq. A3 for i = 1,2 
we find dependences aoi{H), and substitute them into Eq. A4. In accord with 
the discussion at the beginning of the Section III, after Eq. (42), the Lagrange 
multiplier A2 in Eq. A3 plays a role of the " resting" tension in the layer 2 due to 
a finite adhesion energy Ead of the layer to the glassy wall of the pipette. 



2. Semi-equilibrium conformation, P ^ 

In order to find the minimum of (56) we consider NQi, agj and H as the five 
independent variables. Hence, we minimize F with respect to all of them, and 
thus obtain the following equilibrium equations: 

dF/dN^ = : /i - /o - (T/2)(ai - aoi) - {T/2)ao sin 9 = A3 - Aiai; (A5a) 
dF/dN2 = : f2-fo + Ead- {T/2) {a^ - 002) - (T/2)ao sin 6 = 

-A3 - A2a2 ; (A5b) 

dF/dai = : dfi/dai - T/2 = -\i ,i = 1,2; (A5c) 

dF/dH = : Y [9fi/dH - d{T/2)/dH{ai - aoi)] - 

i=l,2 

ao{N2 - No2)d{Tsme)/dH = (Ai + \2)dS{H)/dH . (A5d) 

The system of equations Eqs. A5a-A5d is solved numerically under condition Eq. 
57, after making a transformation to the dimensionless variables Eq. 35 and after 
normalization of the numbers Ni by A^o in accord with Eq. 54. 

3. Equilibrium (relaxed) conformation, P / 

Equations minimizing the free energy F in Eq. 59 take the form: 
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dF/dNi 


= : 


fi — fo — — AiOi ; 


(A6a) 


dF/dN2 


= : 


/2 - /o + Ead - T{a2 - ao2) - Tqq sin © = -X2a2 ; 


(A6b) 


dF/dai 


= : 


dfi/dai = -Ai ; 


(A6c) 


dF/da2 


= : 


df2/da2-T=-X2; 


(A6d) 


dF/dH 


= : 


NM/dH - N2{dT/dH){a2 - aos) - 





1=1,2 



ao{N2 - No2)d{T sin 9) /dH - (Ai + X2)dS{H) /dH . (A6e) 

The system of equations Eqs. A6a-A6e is solved numerically after making a 
transformation to the dimensionless variables Eq. 35 and Eq. 54. 

APPENDIX B: ELASTIC MODEL OF MSCL CHANNEL 

1. Phenomenological theory 

We consider a simple phenomenological model of MscL channel (Sukharev, 
Markin, 2001). The channel is presented as elastic cylinder; it's stiffness is differ- 
ent in closed and open conformations. The free energy of the closed channel, 
is written as: 

Fc = ^{A - Acf - T{A - Ac) , (Bl) 

here Be is elastic modulus of the channel in the closed conformation; T - tension 
transmitted through the lipid bilayer; Ac - is the cross-section area of the closed 
channel in the absence of tension in lipid bilayer; A - is the channel cross-section 
area at a given tension T. Transition to the open conformation of the channel 
is associated with stretching the channel area, rearrangement of transmembrane 
domains (helices) (Sukharev et al., 2001), (GuUingsrud and Schulten, 2003) and 
requires additional work. The work of opening the channel is done by tension in 
hpid bilayer. The free energy of the channel in the open state, Fq, equals: 

Fa = ^{A - Aof + Eo- T(A - Ac) , (B2) 

here Bq - elastic modulus of the open channel {Bq ^ Be)', Aq - the area of open 
channel at zero tension [Aq > Ac)] Eq - the (free) energy difference between 
open and closed conformations of the channel at zero tension in the bilayer. The 
finite energy Eq reflects the fact that the probabilities to find the channel in 
the open or closed states (Po oi Pc) are not equal (at zero tension the MscL is 
closed). In general, the value of Eq depends on many factors, e.g. the type of 
lipids in the bilayer, thermodynamic conditions of the surrounding solution, and 
also on the hydrophobic mismatch energy, F^p (Hamill and Martinac, 2001). The 
latter is different for open and closed conformations of the channel, because MscL 
opening is associated with decrease of it's hydrophobic thickness (S. Sukharev, 
Maryland University, personal communication, 2003). For the given (open or 
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closed) state of the MscL channel, Flp contributes to the total free energy of 
this conformation [Fq or Fq) and thus defines the value of Eq- At tension T 
the minima of the parabolas Fc and Fq define the free energies of the closed and 
open states, and the corresponding areas, Aq or, Alq are the areas occupied by 
closed or open channel in equilibrium: 

Al = ^ + Ac, Al = ^ + Ao, (B3) 

Fc{Al) = , FoiAl) = - TiAo - Ac) + Eo . (B4) 

The equilibrium distribution of open/closed channels, no /t^c i at a given tension 
in membrane depends on the energy difference between open and closed confor- 
mations of the channel, Fq — Fc : 

^ = ^ - exp{-(Fo - Fc)/kBTo} , (B5) 
nc F'c 

here no and nc are fractions of open and closed channels, respectively; ks- Boltz- 
mann constant; T - temperature. The intersection point of the energy parabolas 
Fc{A) and Fo{A) is the barrier energy for transition between closed and open 
states, 

Fc{A,) = Fo{A,) ^ F,,M,) . (B6) 

By comparing (Bl) and (B2), one finds that the area {Ac < < Aq) does 
not depend on the tension, T. The activation energies for closing and opening of 
the channel, Eact,c and Eact,o- 

Eact,c = FbariA,) ~ Fo{Aq) , Eact,o = -^bar(^6) " -^c(^c) i^"^) 

define the kinetics of transitions between closed and open states. The activation 
barriers for closing and opening of the channel equal: 

Eact,c = F,ar{Ab) - Fo{AI) ^ ^{Aq - A,f + T{Ao -A,) + ^- (B8) 

EactAT) = - Acf - T{A, -Ac) + ^. (B9) 

Solving Eq. (B6) with respect to A^ and substituting the result into Eq. (B8), 
one finds dependence £'oct,c on Eq, shown in Fig. 5 for the different values of 
the tension T. As could be discerned already in Fig. 4, Eact,c decreases with the 
increase of Eq. 

The equilibrium distribution of closed and open channels (B5) is reached 
within the characteristic time Teqv' 

re,v = {ko + kc)-\ (BIO) 

here ko and kc are the rates of opening and closing of the channels. The value 
of the rate constant for the open to closed transition, kc , is the key factor for 
the channel closing induced by inter layer slide (as — > 1 ): 

kc = kexp{-Eact,c/kBTo} ■ (BU) 
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2. The energy of hydrophobic mismatch 

The mismatch of hydrophobic thickness of the protein and that of the sur- 
rounding Upid bilayer makes hpids adjust to the protein length. This induces local 
deformation of lipid molecules. The deformation profile, d{r) , can be presented 
in a simple exponential form (Ben-Shaul, 1995): 

d{r) = do + {dp - do) exp {-(r - rp)/ro} , (B12) 

here: r -distance from the protein center in the membrane plane; -radius of 
the protein (protein is modeled as a cylinder); rfg - the (equilibrium) thickness 
of the bilayer hydrophobic core in the absence of protein; dp - is the thickness of 
protein's hydrophobic region; vq - the characteristic scale of the deformation. For 
the sake of simplicity let us consider a flat bilayer. The energy (per unit area) 
of elastic deformation of the bilayer by a cylindrical protein, f^p, to the lowest 
order (neglecting tilt of lipid chains) can be written as (Fournier, 1999): 

hp^^id-dor + ^iVdf, (B13) 

here is the dilation modulus, ~ 2ir^/(2Zo)^; Kg - is the modulus, charac- 
terizing the energy cost of producing a gradient of bilayer thickness (it includes 
the energy of increasing the area of chain- water interface (Fournier, 1999)). 
Wc consider the concentration of the proteins to be small and calculate the hy- 
drophobic mismatch energy, Fip, per channel as follows: 

/•oo 

Flp^ / fLp27rrdr. (B14) 

The result of the integration is: 

Flp = \{Karl + K^) + l) {dp - d^f . (B15) 

Assuming i^^ ~ 15 • IQi^^erg / cm^ (Hamill and Martinac, 2001), Kg ~ 35 ergjcrv? 
(Fournier, 1999), Tq ^ lOA (Ben-Shaul, 1995), Vp ^ 25 A (Sukharev et al, 1999), 
we obtain the following estimates for Flp for two different values of hydrophobic 

mismatch, dp — do : 

1) for \dp — do\ ~ 2A, Flp ~ 2.3A;bTo at room temperature Tq; 

2) for \dp - do\ ~ 6A, Flp ~ 21 kpTo; 

here ks is the Boltzmann constant, Tq is temperature. The value of dp — do 
depends on the type of lipids (i.e. the length of hydrocarbon chains) and on the 
state of the protein channel. 
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TABLES 



TABLE L Numerical results: initial equilibrium at P = 



Ead/et = 0.03 


Ead/et = 0.06 


Ead/et = 0.1 


Ead/et = 0.3 


Z2/Z0 = 1.0077 


Z2/ZQ = 1.0158 


Z2/ZQ = 1.027 


2:2/^0 = 1-093 


n2 = 0.9923 


n2 = 0.9844 


712 = 0.9737 


n2 = 0.9151 


F/{Noet) = 0.03 


F/{NQet) = 0.0596 


F/{Noet) = 0.0988 


F/{Noet) = 0.2878 


z-i_ = ZQ = 0.8291, 


m = 1.000 


for all values of 


Ead/et- 
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Figure Legends 



Fig. 1. Sketch of a membrane in the experimental patch-clamp settings. The 
numbering 1,2 of the monolayers within a bilayer membrane is indicated. 

Fig. 2a. Normalized free energy F/{etNo) in the semi-equilibrium Fdisc and 
complete equilibrium Fconn state versus dimensionless pressure difference p (see 
Eq. 63 ); e./e* = 0.6. 

Fig. 2b. Normalized numbers of lipid molecules in the curved parts of the 
membrane monolayers ni=i,2 in the semi-equilibrium nf^"^ and complete equilib- 
rium 77,?°"" state (see Eq. 54 ) versus dimensionless pressure difference p. Other 
parameters are as in Fig. 2a. 

Fig. 2 c. Normalized curvature of the bilayer membrane inside the pipette 
in the semi-equilibrium hdisc/ho and complete equilibrium hconn/ho state versus 
dimensionless pressure difference p; h/ho = R/r, where R and r are the pipette 
and the curvature radii respectively. Other parameters are as in Fig. 2a. 

Fig. 3 a. Normalized areas per lipid molecule at the external surfaces of 
the monolayers inside the pipette in the semi-equilibrium zf'^'^/zQ and complete 
equilibrium zf^'^ / Zq state versus dimensionless pressure difference p. Other pa- 
rameters are as in Fig. 2a. 

Fig. 36. Dots: Normalized lateral tension Ti in the monolayer 1 in the semi- 
equilibrium state versus dimensionless pressure difference p. Solid line: Normal- 
ized curvature hdisd ho of the membrane inside the pipette in the semi-equilibrium 
state is taken from Fig. 2 c for a convenience of comparison. Other parameters 
are as in Fig. 2a. 

Fig. 4. Sketch of the free energy of MscL as function of the pore area A in 
the elastic model (Sukharev and Markin, 2001). Sohd hues - zero tension T — 0. 
The grey lines indicate free energy at finite tension T. The less intensive are the 
lines - the greater is the tension. 

Fig. 5. The activation energy (barrier) for channel closing, Eact,c as a function 
of the free energy difference between open and closed channel conformations (at 
zero tension) Eq. Solid black line corresponds to zero tension T = 0; the gray 
lines correspond to the tensions T = Qdyn/cm and T = 12dyn/cm (less intensive 
line is for the greater tension). The energies are expressed in units of ksTo. 
Other parameters: Be = 0.22A;bTo; Bq = 2.2A;bTo; Ac = 13nm^; Aq = SOnm'^ 
(Sukharev, Markin, 2001). 
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FIG. 5. 
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